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We study the zitterbewegung of a heavy hole in presence of both cubic Rashba and cubic Dres- 
selhaus spin-orbit interactions. On contrary to the electronic case, zitterbewegung does not van¬ 
ish for equal strength of Rashba and Dresselhaus spin-orbit interaction. This non-vanishing of 
zitterbewegung is associated with the Berry phase. Due to the presence of the spin-orbit coupling 
the spin associated with the heavy hole precesses about an effective magnetic field. This spin pre¬ 
cession produces a transverse spin-orbit force which also generates an electric voltage associated 
with zitterbewegung. We have estimated the magnitude of this voltage for a possible experimental 
detection of zitterbewegung. 
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I. INTRODUCTION 

According to Schrodinger[l], the interference be¬ 
tween two branches of a free Dirac spectrum induces 
a quivering quantum motion, usually known as zit¬ 
terbewegung (ZB). In principle, ZB is not a pure rel¬ 
ativistic phenomenon. In recent years, it has been 
shown that ZB could exist in a plethora of non- 
relativistic physical systems [2] including narrow-gap 
semiconductors [3], spin-orbit coupled low-dimensional 
systems[4-10], graphene[ll 15], carbon nanotube[16], 
topological insulator[17], superconductors[18], sonic 
crystal[19], photonic crystal[20], optical superlattice[21], 
Bose-Einstein condensates [22-25], dichalcogenide mate¬ 
rials like M0S2P6, 27] etc. 

The existence of spin-orbit interaction (SOI) [28] 
in systems like two-dimensional electron/hole gas 
(2DEG/2DHG) is an interesting topic of contemporary 
research due to the promising spintronic applications [29- 
33]. In addition, these systems exhibit various fun¬ 
damental physical phenomena such as zero field spin- 
splitting[34], spin Hall effect [35-39], persistent spin 
helix[40, 41] etc. The absence of structural and bulk 
inversion symmetries in semiconductor heterostructures 
induce Rashba[42, 43] and Dresselhaus [44] SOIs, respec¬ 
tively. The functional dependence of both SOIs on mo¬ 
mentum in 2DEG and 2DHG are different. In the case 
of a 2DEG, they are linear in momentum whereas for the 
2DHG both SOIs are cubic in momentum. Systems with 
cubic SOIs are particularly important as large spin Hall 
conductivity can be achieved in those systems [45-47]. 
Apart from p-doped semiconductor heterostructures, the 
existence of cubic SOIs have been confirmed experimen¬ 
tally in a 2DEG at the surface of SrTiOa[48] and in a 
2DHG in a strained-Ge/SiGe quantum well[49]. 

The problem of ZB of an electronic wave packet in a 
2DEG including both Rashba and Dresselhaus SOIs is 
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well studied. To the best of our knowledge, a less effort 
has been devoted in searching ZB in a spin-orbit coupled 
2DHG. In this work, we study the ZB of the center of 
a Gaussian wave packet in a 2DHG in presence of both 
cubic Rashba and Dresselhaus SOIs. We find that the 
ZB exhibits transient behavior. The amplitude of ZB is 
related with the Berry connection. Unlike the case of a 
2DEG, the ZB survives for equal strength of the Rashba 
and Dresselhaus SOI. We relate this non-vanishing of ZB 
with the behavior of the associated Berry phase. The 
spin associated with the heavy hole precesses about a SOI 
induced effective magnetic field. A transverse component 
of the spin-orbit force is generated as a result of this 
spin precession. Using this spin-orbit force picture, we 
have also estimated the magnitude of an induced electric 
voltage associated with ZB. 

The rest of the paper is organized as follows. In section 
II, we briefly discuss the preliminary informations about 
the physical system. In section III, we consider the tem¬ 
poral evolution of an initial Gaussian wave packet and 
discuss various features of the ZB in position and veloc¬ 
ity. The time evolution of the spin and the spin-orbit 
force picture for ZB are presented in section IV. We pro¬ 
vide asymptotic expressions for ZB and possible detec¬ 
tion scheme in section V. The main results are summa¬ 
rized in section VI. 


II. PHYSICAL SYSTEM 

Let us start with a brief description of the model 
Hamiltonian associated with the spin-orbit coupled 
2DHG. The complicated dynamics of holes at the top 
most valence band of a III-V semiconductor is character¬ 
ized by the 4x4 Luttinger Hamiltonian [50]. However, 
strong confinement in a p-doped III-V quantum well es¬ 
sentially leads to a large splitting between the heavy hole 
(HH) state (|3/2,±3/2)) and the light hole (LH) state 
(|3/2,±1/2)). At low temperature, it is assumed that 
only the HH states are occupied when the density is low 
enough. Now, one can proceed without considering the 
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LH states since a significant contribution to the trans¬ 
port properties near the Fermi energy comes from the 
HH states. Hence, it is possible to obtain an effective 
2x2 Rashba Hamiltonian by projecting the 4x4 Lut- 
tinger Hamiltonian onto the HH states [45]. Additionally 
the Dresselhaus SOI originates due to the bulk inversion 
asymmetry of the host crystal[51]. 

Now, the single particle dynamics of a HH is governed 
by the following Hamiltonian 
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where a(k, 0) = d/(y/n)e~^ d ( k ” k °) with d (ko) is the 
width (wave vector) of the initial wave packet, gu = 
i e -i{2e-<i> k ) anc j Wk _ ( e + _ e ~)/2 h. However, the wave 
packet was initially polarized along +z direction but it 
acquires a component along —3 direction as time goes on. 

By taking the inverse Fourier transformation of Eq. (5), 
we obtain the following expression for the wave packet in 
the momentum space at time t 
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A. Calculation of the expectation values 


where a± = a x ± ia y with cq’s as the usual Pauli spin 
matrices and p± = p x ± ip y with p is the momentum, 
m* is the effective mass of the heavy hole and a (/3) is 
the strength of the Rashba (Dresselhaus) SOI. Note that 
a can be tuned by an external gate voltage but fj is a 
fixed material dependent quantity. Here, Pauli matrices 
represent an effective pseudo-spin with spin projection 
±3/2 along the growth direction of the quantum well. 

The eigenvalues and eigenstates of the Hamiltonian are 
respectively given by [52], 

^=S? +AA * fc3 ’ (2) 

and 

where A = ±, A g = ^ a 2 ± /3 2 — 2a/3 sin 29 with 9 = 
ta,n~ 1 (k y /k x ) and <f>k = tan _1 [(afc x — /3k y )/(ak y — /3k x )]. 
Note that the energy spectrum (Eq. (2)) is highly 
anisotropic described by the factor Ag. The spin split¬ 
ting between the HH branches at a given wave vector can 
be obtained as ej/ - e^=2Agk 3 . 

III. TIME EVOLUTION 

Here, we seek to see the time evolution of an initial 
wave packet representing a HH in the presence of the 
SOIs. Now, applying the time evolution operator U(t) = 
e ~iHt/h on thg initial Gaussian wave packet polarized 
along ~t~z axis 

4, ( r ’°) = 2^ / d2 MM)e ?k ' r ( o 'j , (4) 

we find the wave packet at a later time t as 

tf(r ,t) = ^~ I d 2 k a(k, 0)e ik r e~ i ^ rt 
2ir J 

x | cos(wkt) o ) + sin(a;kt) ^ ^ |,(5) 


Here, we would like to calculate the expectation val¬ 
ues of position and velocity operators explicitly. The 
expectation value of any physical operator O is defined 
via ( O } = / d 2 fc$^(k, t)0$(k, t). After a lengthy calcu¬ 
lation, one can obtain the average value of the position 
operator as 

<r(*)> = <r( 0 )) + ^® t _ 1 f d 2 fc|a(k, 0 )| 2 

m* 2 J 

x (v k (20-</ fe )){l-cos(2 Wk t)}, (7) 

where (r(0)) and (k(0)) can be determined from the ini¬ 
tial condition. The last oscillatory term in Eq. (7) repre¬ 
sents the phenomenon zitterbewegung , the frequency of 
which is governed by the spin-split energy differences i.e 

w k = (4 - e k)/( 2 ^)- 

Now, setting (r(0)} = 0 and (k(0)) = k 0 = koy, we ob¬ 
tain the following expressions for the expectation values 
of the position operator 

rj r2n roo 
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where q = kd , a o = kod , Q = y/\ ± r/ 2 — 2rj sin(20) with 
r] = /3/a, and = 2 aq 3 Q/(hd 3 ). Note that when fj = 0, 
there will be no oscillatory term in Eq. (9) since the in¬ 
tegral J 2k e 2a ° qsin 8 cos 9d9 is exactly zero. However, for 
a finite /3, the angular integration can not be done an¬ 
alytically. Nevertheless, one can show numerically that 
the second term in Eq. (9) is negligibly small compared 
to the first term. In this way, one can say that the ZB 
essentially occurs in a direction perpendicular to the ini¬ 
tial wave vector. Henceforth, wc will consider only the x 
component of the observables. 
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The x component of the velocity operator can be ob¬ 
tained using the Heisenberg equation v x = [x, H]/(ih) 
as 

Px ( 2 2 \ 

"• = ■^‘ T « + WV+ p -- ,T - p +) 

- +P-) + ^-( 2 P 2 +P+)}- ( 10 ) 

After a straightforward calculation, we find the expec¬ 
tation value of v x as 

p 27 T poo 
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where v r = a/(hd 2 ) and g{0) = g 2 + 3 — 4?ysin(20). 



B. Non-vanishing of zitterbewegung for a = 0 

A close inspection of Eqs. ( 8 ) and (11) reveals that 
the ZB appearing in either position or velocity does not 
vanish for a = 0- This result strongly contradicts the 
electronic case in which it was argued [ 6 ] that ZB would 
vanish for equal strength of Rashba and Dressclhaus SOI 
because of the existence of an additional conserved quan¬ 
tity <j x — <J y . But for the case of a 2DHG no such con¬ 
served quantity exists. However, here we try to relate 
the non-vanishing of ZB with the Berry phase associated 
with the energy spectrum. The Berry connection is given 
by Ak = (V'ksN^IV'ks); where ip£ s denotes the spinor 
part of the wave functions given in Eq. (3). An explicit 
calculation will yield 

1 ( cJ 2 — B 2 \ 

A k = — ^2+ P j(k x y-k v x). (12) 

The corresponding Berry phase can be obtained by 
taking the line integral of the Berry connection as 7 = 
f A k • dk. Using the residue theorem of complex inte¬ 
gration we find the Berry phase as 

ct 2 — 0 2 

7 “ 2 ’' + ’ r i^vyy < 13 > 

Eq. (13) clearly shows that the Berry phase does not 
vanish for the case of a = B- It has been found[53] for 
the 2DEG case that the Berry phase becomes zero when 
the strengths of Rashba and Dresselhaus SOI are equal. 

C. Transient zitterbewegung in position and 
velocity 

In Fig. 1 we have shown the time variation of the 
square of the autocorrelation function (|A(f)| 2 ), expec¬ 
tation values of position ( x{t)) and velocity {v x {t)) op¬ 
erators for a fixed ao = 10 and different g = 0, 0.5, 


FIG. 1: Plots of (a) square of the autocorrelation function, 
(b) position, and (c) velocity expectation values for ao = 10. 
Solid, dashed, and dotted line represent g = 0, 0.5, and 1, 
respectively. Here, t is considered in units of 10 -3 ftd 3 /(2a). 

and 1. The autocorrelation function mainly measures the 
overlap of the time evolved wave packet with its initial 
counterpart and is defined as A(t) = (^(r,f)|\H(r,0)) = 
f d 2 fc $i(k, t)<h(k, 0). Now using Eq. (6), it is straight¬ 
forward to find A(t) in the following form 

A(t) = i e~ a ° j d9dqqe- q2+2aoqsin8 cos 

(14) 

Note that for a finite a / 0, / 0, the integrals over 

q and 6 appearing in Eqs. (8), (11), and (14) can not be 
done exactly and hence numerical treatment has to be im¬ 
plemented to see the explicit time dependence of |A(f)| 2 , 
\x(t)}, and (v x (t)). From Fig. 1(b) and Fig. 1(c) it is 
evident that the ZB, appearing in position and velocity, 
shows transient oscillations. The time scale associated 
with ZB can not be extracted analytically. However, the 
transient character of ( x{t )) and (v x (t )) may be explained 
from the behavior of the autocorrelation function which 
is shown in Fig. 1(a). |A(t)| 2 exhibits damped oscilla¬ 
tory character. Fig. 1 clearly shows that the oscillations 
associated with ZB stop when the autocorrelation func¬ 
tion nearly dies out. Now the effect of an enhancement 
of B on ZB is twofold. First, it introduces a phase in the 
oscillations and second, the amplitude of ZB decreases 
with the increase of /? as evident from Fig. 1(b). 

IV. SPIN DYNAMICS AND SPIN-ORBIT 
FORCE 

Here, we are interested to study the time dependence 
of the expectation values of the effective spin operator 
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FIG. 3: Plots of the expectation values of spin-orbit force op¬ 
erators for ao = 10. Solid, dashed, and dotted line represent 
77 = 0, 0.5, and 1, respectively. Here, t is considered in units 
of 10~ 3 hd 3 / (2a). 


FIG. 2: Plots of the expectation values of the spin operators 
for do = 10. Solid, dashed, and dotted line represent r/ = 
0, 0.5, and 1, respectively. Here, t is considered in units of 
10~ 3 hd 3 /(2a). 


and the spin-orbit force. Due to the existence of SOIs, 
the spin S = (3/2 )hcr associated with the HH precesses 
about an effective magnetic field Cl in accordance with 
the following equation 

$ = Cl x S. (15) 

dt 

The components of Cl can be obtained as 
= ^[ap y {3pl -p 2 y ) - 0p 2 p x ] 

and 

2 

tty = ¥ lap x (3pl - pl) - Pp 2 Py\- 
The general solution of Eq. (15) is given by 

S = (n ■ So )fi + [So — (n • So)h] cos(f It) 

+ (n x So) sin(fli), (16) 


where Cl = \JCl^. + Cl 2 , So is determined from the initial 

condition, and the unit vector n is defined as n = Cl/Cl. 
Without any loss of generality, we choose the initial con¬ 
dition such that h • So = 0 or equivalently, So = S z z 
since n lies in the xy plane. So, one can obtain the com¬ 
ponents of the effective spin operator at a later time t as 
S x (t) = (Cl y /Cl)S z sm(Clt), S y (t) = — (Cl x /Cl)S z sin(fii), 
and S z (t) = S z cos(Clt). We find the expectation values 
of the spin components at time t as 
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x cos(Cl q t). 
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(18) 


Let us now derive the expression of the spin-orbit 
force. The spin-orbit force is defined as F = m* r, 
where r can be evaluated from Heisenberg equation r = 
— [[r , H], H]/h 2 . Now, it is straightforward to obtain F 
in the following operator form 


F = 


2 m*p 

~h7~ 


7= 2 r 


(3a 2 + P 2 )p 2 - 8a(3p x p y (p x z)a z (l9) 


Using Eq. (6) we calculate the expectation values for the 
the components of the spin force operator as 
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2m* a 2 2 f 2 * 

Trh 2 d 5 J 0 
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( 20 ) 


In Fig. 2 and Fig. 3, we have shown the time de¬ 
pendence of the expectation values of spin and spin-orbit 
force operators. From Fig. 2 it is evident that ( S x (t )) 
vanishes in the absence of the Dresselhaus SOI. When 
0 = 0 the angular integrations in Eqs. (17) and (18) can 
be done exactly which gives a vanishing result for ( S x ). 
The other components of the spin operator i.e. ( S y (t )) 
and ( S z (t )) can be found in the following forms 

p OO cy . 

(S y (t)) = 2>he~ a ° j dq qe~ q2 1 3 (2a 0 q) sin (— q 3 t)(21) 

























5 


and 

2 2 / 2a n \ 

(S z {t )} =3 he~ a o J dq qe~ q / 0 ( 2 a 0 g) cos (^3 q 1 2 3 tj ( 22 ) 

where I n (z) is the modified Bessel function of first kind. 
Eqs. (21) and (22) indicate that the spin precesses in the 
yz plane. For f3 = 0, one can find that ( F v (t )) = 0. In 
this case, ( F x [t)) takes the following form 

(F x (t)) = F 0 e~ a ° J dq q 6 e~ q ~I\ (2aoq) cos {^^q 3 tp. 3) 

where Fq = 12 m*a 2 /(h 2 d 5 ). Analyzing Eqs. (21)-(23), 
one can conclude that the spin precession in the yz plane 
essentially generates a transverse spin-orbit force along x 
direction. This fact is depicted in Fig. 3. 

The scenario changes significantly in the presence of a 
finite /?. In this case the angular integrations in Eqs. (17)- 
(18) and Eq. (20) can not be done exactly. A numerical 
calculation shows that ( S x (t )) survives for j3 ^ 0 but 
its magnitude is small compared to that of ( S y (t)) and 
(. S z (t )). Thus the spin starts to precess in the entire space 
as one switches on /?. As a result, a finite y component 
of the spin-orbit force is generated. But the amplitude 
of (F y ) is negligibly small in comparison to that of (F x ) 
as evident from Fig. 3. So, one can safely consider that 
the spin-orbit force acts along the x direction provided 
the initial wave packet is injected along the y direction. 
In recent past [5], the origin of ZB was explained in the 
light of the spin-orbit force and it was argued that ZB 
occurs in a direction along which the spin-orbit force acts 
on the wave packet. By observing Eqs. ( 8 ), (11), and 
(23), we conclude that our results are consistent with 
that prediction. 

V. PERSISTENT ZITTERBEWEGUNG AND 
POSSIBLE DETECTION 

The ZB, manifested in the expectation values of 
the physical observables undergo transient behavior. 
However, the persistent behavior of ZB can be ob¬ 
tained in the limit ao > 1 i. e. the initial wave 


packet is sharply peaked in the momentum space. 
For /3 = 0, one can obtain the following approxi¬ 
mate expressions for the expectation values: (x) — 
(3/2fc 0 )[l — cos(w 0 t)], (v x ) ~ (3afcg/fi)sin(wot), (S y ) ~ 
(3h/2) sin(wot), (S z ) ~ (3H/2) cos(wot), and (. F x ) ~ 
/ocos(wot) where wo = 2ako/h and /o = 6m*a 2 kl/h 2 . 
Due to the presence of SOI, an effective Lorentz-like 
magnetic field [54] can be obtained in the following form 
B s = 2m* 2 a 2 p' L g(0)a z z/(eh 7 ). In the ao > 1 limit and 
when /3 = 0 we obtain (B s ) ~ B 0 cos(wot) where B 0 = 
6m* 2 a 2 k^/(eh 3 ). Interestingly, we find (. F X )/(B S ) = 
ehko/m* i.e. F x plays a role of the Lorentz force in the 
ao 1 limit. Now, the time dependent magnetic field 
B s would generate an electric voltage of the following 
form V = Vo sin(woi) where Vo = woABq with A as the 
area of the sample. For typical material parameters a 
GaAs quantum well i.e. m* = 0.45ro e , a = 10 -28 eVm 3 , 
A = 10 4 nm 2 , and ko = 10 8 m _1 we obtain ojo — 3 x 10 11 
rad/s, B 0 ~ 0.14 T, and Vo ~ 4.4 x 10 -4 V. In principle, 
one can measure Vo to confirm the signature of ZB exper¬ 
imentally. As introduced recently[10] this kind of electric 
voltage associated with ZB can be also obtained in the 
case of a 2DEG with time dependent linear Rashba SOI. 


VI. SUMMARY 

In summary, we have investigated the zitterbewegung 
of a wave packet in a 2DHG in the presence of both 
Rashba and Dresselhaus SOIs. The zitterbewegung man¬ 
ifested in position, velocity, and spin appear as transient 
oscillations. For equal strength of Rashba and Dressel¬ 
haus SOIs we find zitterbewegung does not vanish, a 
direct contradiction to the electronic case. This non¬ 
vanishing of zitterbewegung has been interpreted as a 
consequence of the Berry phase. The pseudo-spin associ¬ 
ated with the heavy hole precesses about a SOI induced 
effective magnetic field. As a consequence of this spin 
precession, a transverse component of the spin-orbit force 
appears which in turn induces zitterbewegung in posi¬ 
tion and velocity. The transverse spin-orbit force also 
generates an electric voltage. We have concluded that an 
indirect signature of zitterbewegung can be unveiled by 
measuring this voltage. 
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